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ABSTRACT 

Preston-tube measurements obtained on the Arnold Engineering Development 
Center (AEDC) Transition Cone have been correlated with theoretical shin 
friction coefficient in transitional and turbulent flow. This has been done 
for the NASA Ames U-Ft Transonic wind Tunnel (11 TWT) and flight tests. 
During analyses of the data,, errors were discovered in the Preston-tube 
measurements. The errors in the 11 TWT data are -minor and were easily 
corrected. The errors in the flight data, however, are muchjarggr and 
random. The source(s) of these errors are unknown. We suspect they are due, 
at least in part, to twisting of the probe during a traverse. A procedure has 

been developed to correct for these errors and has been successfully applied 
to the flight data. 

The developed semi -empirical correlations of Preston-tube data have been, 
used to derive a calibration procedure for the 11 TWT flow quality. This 
procedure has been applied to the corrected laminar data, and an effective 
freestream unit Reynolds number is defined by requiring a matching of the 
average Preston-tube pressure in flight and in the tunnel. In contrast to 
previous studies of the effect of tunnel noise on boundary layer transition, 
tbds study finds that the operating Reynolds number is below the effective 
value required for a match in laminar Preston-tube data. The distribution of 
this effective Reynolds number with Mach number correlates well with 
freestream noise level in this tunnel. Analyses of transitional and turbulent 
data, however, did not result in effective Reynolds numbers that can be 
correlated with background noise. This is a result of the fact that vorticity 
fluctuations present in transitional and turbulent boundary layers dominate 
Preston-tube pressure fluctuations and. therefore, mask the tunnel noise 
effects. So, in order to calibrate the effects of noise on transonic wind 
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tunnel tests only laminar data should be used, preferably at flow conditions 
similar to those in flight tests. To calibrate the effects of transonic wind- 
tunnel noise on drag measurements, however, the Preston-tube data must be 
supplemented with direct measurements of skin friction. Such data could be 
used in the subject procedure to define an equivalent flight unit Reynolds 
number for a given tunnel setting jadJich would result in a matching of the 
average, laminar skin friction_ 
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CHAPTER I 

DEFINITION OF THE PROBLEM 

Since the transonic wind tunnel became operational at NASA Langley In the 
late l&40's, there has been a need for a procedure to calibrate the effects of 
wall -generated noise on the tunnel flow quality. As noted by Dougherty and 
Stelnle (1), the primary Indicators of flow quality in a wind tunnel are 
variations of: (1) Mach number, (2) flow angularity within the empty test 

section, and (3) the Reynolds number at which transition from -laminar to 
turbulent -flow occurs on models. Variations In Mach number and flow 
angularity can be calibrated with, conventional Pitot-static probes and yaw 
meters, e.g. see Reed, et al. (2) And In the case of low-speed wind tunnels, 
the Reynolds number at which the drag coefficient of a sphere equals 0.30 can 
be used to define a turbulence factor_(TF), as described by Pope and Harper 

(3). An effective" unit Reynolds number for a given low-speed tunnel can 
then be defined by 


Re m.eff * < TF > R <*m 

However, when Mach number exceeds about 0.35. compressibility, effects 
cause the classical turbulence factor to become increasingly erroneous and 
therefore not useful. Recently, Miller and Bailey (4) have reviewed the 
status of knowledge concerning the drag of a sphere at transonic speeds. Even 
today, the precise variation of sphere drag with Mach number and Reynolds 
number is not well defined. Thus, the classical turbulence-sphere method is 
not applicable to the calibration of transonic wind tunnels. 
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In recognition of the non-appl icabl Itty of a sphere for defining a 
turbulence factor for tests in transonic wind tunnels, NASA, as part of the C- 
5A wind tunnel correlation program (Treon, et al.(B)) employed what is now 
known as the AEDC 10° Transition Cone as a means of defining an adjustment to 
Reynolds numbers on a tunneL-to-tunnel basis. The cone was initially 

developed in the mid-sixties by engineers at Arnold Engineering Development 
Center (AEDC). It has a traversing Pitot probe resting on its surface to 
directly detect boundary-layer transition.! The cone geometry h.as the 
advantage that no shock is generated along the surface at transonic speeds, 
and theceby avoids shock/boundary-layer interactions such as occur on 
airfoils, wings and blunt-nose bodies. A schematic of this cone and some of 
the associated instrumentation are shown in Fig. l. since the cone was 

designed to calibrate the effects of tunnel noise on boundary-layer 
transition, it also has two miniature microphones imbedded in the surface at 
45.7 and 66 cm aft of the nose for noise measurements. Additional description 
of this cone can be found in the papers by Dougherty and Steinle (1) and 
Dougherty and Fisher (7). 

The need for such a calibration device was emphasized when discrepancies 
between numerous transonic wind-tunnel tests of models at ostensibly identical 
flow conditions were observed. A particularly we 1 1 -documented study of 
differences in static aerodynamic data has been obtained with the same model 
of a Lockheed C-5A transport aircraft in three major transonic wind tunnels. 
The results have been reported by Treon, et al. (5) The differences between 
the three different sets of wind-tunnel data were reduced by accounting for 


P 1 ] 3 ' C ?u rSe ! i s not a new measurement technique. In fact, the first 

this^technitju^fo^filght^tests?* (6) 1937 deSCr,bes the ° f 
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"relative" Reynolds number effects between facilities. The AEUC Transition 
Cone was used to define the differences in "relative" Reynolds number. 

As observed by Dougherty and Steinle (2): "These results substantiated 
the need for developing a method for predicting these-corrections to Reynolds 
number to improve extrapolation of wind-tunnel test results to full-scale 
conditions, i.e., a "turbulence factor for transonic tunnels." 

Subsequent to the C-5A correlation program, the cone was tested in 
transonic wind tunnels both in the U.S. and in Europe. At the completion of 
the wind tunnel tour, the cone was tested in flight, Dougherty and Fi sher 
(6). Parallel with the planning of the flight program, NASA focused on using 
the data from the cone in flight, in conjunction with the wind tunnel data, to 
develop a means for defining an adjustment (transonic turbulence factor) to 
Reynolds number on a tunnel -to-fl ight basis. 

Such a Reynolds number will calibrate noise effects on the onset of 
transition, so that by increasing, the flight Reynolds number to that value, 
transition will occur at the same location as in the tunnel. However, 
matching of transition onset is not sufficient practical use, since other 
useful- parameters like Preston-tube pressure and skin friction measurements 
are not necessarily matched by that procedure. 

The objective of this work, on the other hand, is to infer skin friction 
along the AEDC cone using the Preston-tube Impact pressure measurements In 
both wind-tunnel and flight tests and, in analogy with the turbulence sphere 
method, define a procedure whereby an "effective" freestream unit Reynolds 
number can be calculated for a given tunnel setting, but this number now 
represents the freestream unit Reynolds number at which the model tested in 
the tunnel will experience the same average, theoretical skin friction as In 
flight, or, equivalently, will give the same measured values of average 
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Preston-tube pressure. 

The importance of this work lies not only In the calibration of wind 
tunnel flow quality, but also in the general and systematic way of relating 
wind-tunnel flow conditions to actual flight. Thus, the prediction of flight 
level drag will be improyed and the results obtained from wind tunnel tests 
can be directly applied towards the desi gn and development of prototypes. 

The basic approach used in this study to achieve the above-mentioned 
objective is as follows: 

1. Preston-tube measurements are correlated with theoretical skin friction 
along the surface of the AEDC Cone_in laminar, transitional and turbulent 
portions of the boundary-layer flow. This is done for the wind tunnel tests 
as well as the flight tests. 

2. With the two sets of correlations {one set of three correlations for the 
wind tunnel tests, and a second set for the flight tests), the skin friction 
coefficient is equated as well as all other variables and parameters, except 
the freestream unit Reynolds number, Re m . The two sets of correlations are 
expected to have different empiuial coefficients since noise and freestream 
turbulence effects, which are not modeled in the theoretical computations, are 
different. This means that substituting wind tunnel data, which includes 
but ex clude^ Re m , into the flight correlation results in a freestream unit 
Reynolds number that is different from the measured one in the tunnel. This 
derived Reynolds number is therefore the noise-free "effective" Reynolds 
number the tunnel should operate at to get the same P p measurements as in 
flight. 

3. Since correlations for the laminar portion of the boundary layer are 
expected to be different from those in the transitional and turbulent 
portions, the resulting Re m>eff 's may be different in general for each 
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portion. Analysis of these results should reveal the best measure of tunnel 
flow quality. 

An earlier report by the authors (8) presented the results obtained- from 
the analysis of laminar data, the— adjustments to freestream unit Reynolds 
numbers, AR eff » were found to range between 3.0 and 6.5%. The findings of 
that report are summarized at the beginning of Chapter IV. The present report 
focuses on the analysis of transitional and turbulent^data. 


CHAPTER II 

SURVEY OF RELATED LITERATURE 

Since Pceston-tube pressures are by definition total pressures near the 
wall, the classical law-o£=ihe=*a11 can be used to relate these pressures with 

wall shear stress. The law-of-the-wall can be expressed in the following 
general form: 


U + - F x (Y + ) . 


(I) 


Usl^g the definitions of U + and Y + , Eqn. (1) can be written as 


where U is the velocity parallel to the wall at the normal distance Y. 
Associated with the Preston-4;ube measurement of total pressure, P p , is a 
velocity. Up, at a height Y^f. In other words, there exists a streamline 
entering the probe face, Y eff units above the wall, at which the theoretical 
total pressure In the undisturbed boundary layer flow equals the total 
pressure measured by the probe, P p , Fig. 2. This "effective" probe center or 
height concept was introduced by Preston (9) in 1953. The corresponding 
theoretical velocity at this height Is denoted by U D . 

r 

Thus, at the effective height, Eqn. (2) Is written as follows 


U P U T Y *ff 

_2_ - r / T eff \ 

TT F l ( — — > 

T W 


Multiplying both sides by U Y „/v gives 

t err w 


( 3 ) 
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U pV.f = Meff F / U T Y eff ) « F / U t Y eff v 

v„ V r l ' "“v > * \ —z )* 


or, alternatively. 


U Y „ 
t eff 


w 


* F* ( 


U P Y eff 


). 


( 4 ) 


Eqn. (4) Is the general form of the correlation between Preston-tube 
measurements and skin friction. 


Incompressible-Flow Correlations 

According to Preston (9), the British engineers Stephens and Haslam (10) 
suggested in 1938 that It should be possible to use the data from a Pitot tube 
traversed along a surface to infer skin friction. Apparently, this Idea was 
not pursued until Preston's work during the early 1950 's. He developed a 
correlation between skin friction and the. total pressure as measured with 
circular Pitot tubes resting on the inside wall of a pipe. In order to 
develop his correlation, Preston used a simplified version of Eqn. (4) by 
making two assumptions: 

1. The flow Is Incompressible and Bernoulli's equation is valid, thus U p 
can be easily related to P p as follows: 


2. The effective center of the circular tube is fixed and coincides with 
Its geometric center, l.e. Y eff * d/2. 

These two assumptions lead to the following relation. 
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r i 


T w d 

4 Vw 


(p - P ) d* 
v p w' 

4 p v 2 
w w 


Using Eqn. (5J as a guide, Preston obtained measurements Inside a pipe flow 
with circular Pitot tubes having four different external diameters but a 
nearly constant ratio of internal to external diameter of 0.6. Pipe Reynolds 
number was varied over the range 10 4 < Re D < 10 5 . Skin friction was 
determined via measurements of pressure drop over a known length of constant- 

diameter pipe, viz., t w = (P 1 - r 2 )d/ 4L. An empirical fit of the data led to 
the following correlation. 

* - 1.396 + 0.875 xj ( 6 ) 

Where the variables are defined as 


(P - P ) f 

»? ■ l°9l0 [-*- —*$ 3, (6a) 

4 p v 

T d 2 

y ! * ’OHIO c 7 — r (6b) 

4 p v 

In 1964, Patel (11) published the results of an extensive set of tests 
with fourteen circular Pitot probes and three different pipe diameters. He 
obtained a more accurate calibration for Preston tubes and established limits 
on the pressure gradient conditions within which his calibration can be used 
with prescribed accuracy. Patel obtained empirical equations for y* = f( x *) 
over three regions of y*. These regions correspond to the ful ly-turbulent , 

the buffer or transition zone, and the viscous-sublayer regions of the 
classical law-of-the-wall. The normal Reynolds number range of Preston-tube 
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measurements in Incompressible flow correspond to the. buffer zone, and for 
this region Patel obtained 

yf = 0.8287 - 0.1381X* + 0.1437(x*) 2 - 0.0060(x*) a , (7) 

where 1.5 < yj < 3.5 or 5.6 < U T d/v w < 55. Patel reported this correlates his 

data to within + 1.5% of t, . 

— w 

In the viscous sublayer region, Patel found his data were correlated by 

y* = 0.5 x* + 0.037, (8) 

when y* < 1-5 or d/v w < 5.6 . In this near-wall region, the classical 
law-of-the-wall exhibits the linear relation 

U + = Y + . (9) 

In order to relate (8), (9), Patel introduced K eff as the normalized effective 
center of a round Pitot tube defined by 

K eff = 2 Y eff /d * ( 10 ) 

Substituting into (9) and using the definitions of x* and y* result in the 
following equation. 


y* a 0.5x| - 0.5 1og 1Q (0.5 (11) 

When this equation is equated with Eqn. (8) and solved for K eff , the result Is 
Keff = 1.3 . 
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The traversing Pitot probes, used during wind-tunnel tests with the AEDC 

transition Cone, are of the flattened or oval-shaped type. Since Patel's 

correlations are for circular Preston tubes, they cannot be applied directly 

to the AEDC Cone tests.. In addition, these tests were conducted at transonic 

speeds, and compressibility effects are expected. With regard to the 

flattened Preston tubes, Quarmby and Das (12) conducted an experimental study 

and calibration of six oval-shaped Preston tubes. When xf > 4.6, they found 

these probes gave exactly the same calibration, relation, between yfand xf as 

was .obtained by Patel (Eqn. 7) if the external height of the probe face is 

used in place of d. At lower values of x£ , the negative displacement of 

effective center caused by wall proximity was larger (= 5%) for the flattened 

probes w<th aspect ratios between 1.5 and 1.9.t The following calibration 

equation correlated the measurements of Quarmby and Das within 1.5% of t . 

w 

y\ = 0.5152 + 0.1693xf + 0.0651 (x*) 2 , 

( 12 ) 

3.38 < x* < 6 . 

The two correlations, Eqns. (11) and (12) make the assumption that the 
effective center of the probe Is fixed.. Preston showed that It is a function. 
of U t d/y w but djd not attem P* to define this function. 

McMillan (13) pursued this point and found for circular tubes that the 
displacement of the effective center is 0.15d (K eff = 1.3) when the probe is 
more than two diameters away from the wall and Is affected by shear flow 
alone. As the probe gets closer to the wall, K ef f decreases. McMillan 
confirmed, therefore, that K e ff is a function of U T h/v w » One can understand 


t This Is consistent with the Idea that flow about the face becomes more 
two-dimensional as aspect ratio Increases and more of the flow passes up 
and over the face rather than around the sides. 


this wall proximity effect by considering that a greater portion of the flow, 
blocked by the probe, will have to lift upward and move over and around the 
probe face as less passes underneath between the probe and the wall. McMillan 
proposed a single curve, Independent of Reynolds number, to correct for wall 
proximity effects on the measured velocity. 

The work done by Patel (11), McMillan (13) and Quarmby and Das (12) leads 

to the conclusion that. In general, K e ff_1s a function of U h/v , Y /h and 

t w y 

w/h (aspect ratio). Since we are talking about a Pitot tube resting on. the 
wall, Yg/h * 0.5. and since the same probe_ls_iiS£d in all the tests, w/h Is 
constant. Therefore, the relation 


U h 

K eff s F n ‘ ( ^ (13) 

w 

seems to describe the actual variation In K eff for Incompressible-flow 
conditions. If this relation is Incoreprated in Eqn. (4), It can be shown 
that can be eliminated, while Eqn. (4) remains in the same form. This 
explains why the assumption of fixed effective probe height has worked well 

for Incompressibl e-flow correlations. 

For compressible-flow correlations, however, Eqn. (13) is expected to be 
different. It will perhaps have the form 

U h 

K eff “ Fn * ( — * MJ. (14) 

w 

In this case,_any attempt to neglect the variation of K eff must show up In a 
greater scatter of data about the developed correlation. 


Compressible-Flow Correlations 

Allen (14) has performed a comprehensive analysis of Preston tubes In 
supersonic boundary layers. He developed a correlation using three 
Independent sets of simultaneous measurements of Preston-tube pressures and 
skin friction via floating-element force balance. These data were obtained 
within flat-plate, turbulent boundary layers and with freestream Mach .numbers 
in the range: 1.6 < < 4.6. Allen selected the same basic parameters as 

Patel; except, he chose to evaluate the fluid properties at a reference 
temperature developed by Sommer and -Short (15), and the velocity U p was 
calculated from.Rp and the_wal 1-pressure P w (=P Q ) using standard compressible 
flow relations.! - 


x 2 * '°9l0 1 


v 


u d 

yf = lo 9l0 < 72- > 


(15a) 


The primes denote properties evaluated at the Sommer and Short reference 
temperature, viz., 


j~ = 0.55 + 0.035 M? + 0.45 ^ . 
e e 


The correlation derived by Allen Is 


(16) 


y* = -0.4723 + 0.74814 x£ + 0.01239 (x£) 2 . 


( 17 ) 


t The details can be found In a report by Allen (16). 


Allen found that the majority of the skin-friction-coefficient data were 
within +1K% to *12% of Eqn. (17). This rather Urge scatter, compared to the 
Incompressible pipe-flow calibrations of Patel and Quarmby and Das, Is at 
least partly associated with the much greater sensitivity and vulnerability of 
floating-element balances to extraneous errors. t 

Obviously, the parameters used by Allen are logical candidates In any 
attempt to correlate the transonic cone data- However, the basic -purpose of a 
reference temperature is to permit use of skin-friction formulas for 
Incompressible flow to estimate compressible skin friction by evaluating fluid 
properties at the reference .temperature. Thus, the resulting reference 
properties represent "average" values across a boundary layer. Whereas, small 
Preston tubes encounter only the flow near the wall. Therefore, It appeared 
to the .authors that properties based simply on the wall temperature would be 


more appropriate. 


Laalnar Preston-Tube Correlations 


A survey of the literature uncovered only one paper, published by 
Prozorov (19) which addresses the. problem of using Preston-tube measurements 
to deduce skin friction in a laminar boundary layer. He obtained surface 
Pitot-probe measurements within low^speed, flat-plate, laminar boundary 
layers. He used several circular and rectangular probes with different aspect 
ratios. Though his data exhibited considerable scatter, he concluded that 
K e ff Is a function of U d/v w for both laminar and turbulent portions of the 


t Allen (17) discussed the various error sources in floating-element 

balances. He has recently suggested an Improved design for this type of 
Instrumentation, Allen (18). 



boundary layer Irrespective of the aspect ratio, which Is Inconsistent wit 



the results of McMillan U3) and Quarmby and Das (12). He also found 

T w d-/p w vf\ (the square of U d/v ) to be a different function of U d/v 
“ " x w p w 

compared to what Preston (9) found. 

His deduction jo£ the laminar correlation is based on a McLaurln series 
expansion _o£_ U p near the wall (since the probe height was small relative to 
the boundary layec thickness) and the conservation equations of mass and 
momentum for__steady, two-dimensional, incompressible flow. The result is the 
following equation. 


T w 


u U . dP 
w p 1 e 

T7T ~ iir 


eff 


eff 


(18) 


Prozorov's correlation takes into account the pressure gradient. The 
theoretical calculations of inviscld static-pressure distribution by Wu and 
Lock (20) for the wind tunnel cases, and the measurements of surface pressures 
in the flight tests show that the pressure gradient in this study is 
negligible. Prozorov claims that his correlation is valid in laminar, 
transitional and turbulent flows provided that the probe is always w1thin_the 
viscous sublayer; a condition which was found to be invalid in this study. 

It can be shown that Prozorov's correlation is equivalent to the 
calibration model used in this study when dP g /dx * 0 and K eff is small. 



Boundary-Layer Transition Computation 

Dhawan and Narasimha (21) developed a method of calculating the 
properties of a boundary layer undergoing transition by preserving the 
essential Intermittency of the flow. Narasimha (22) modified Emmons's (23) 
original function to obtain .an intermittency function described-by 


Y(x) « 1 - e~ A52(x) , a = 0.41 , ( 19) 

where 

x - x R 

S( X ) = — i — “ • (20) 

Here Xg is the transition point defined as the location where the Pitot-tube 
measurements depart from the laminar ones and x Is defined by 

A = V.75 " V.25 * (21 > 

By comparison with numerous other data, including supersonic data, Eqn. (19) 
was shown to be a good approximation to a universal intermittency function for 
boundary layer transition. Therefore, this function is adopted in this study 
to model boundary, layer transition, and its use Is. as will be described in 
detail In Chapter IV. Fig. 3 illustrates a typical y distribution and how It 
changes with X. 




Calibration of Wind-Tunnel Flow Quality 

With the establishment of the fact that freestream disturbances can 
significantly affect transonic wind-tunnel data, an extensive test program was 
begun at the NASA Ames research center In 1971. The AEDC Cone_was tested- in 
twenty-three tunnels between 1971 and 1977. Finally in 1978, it was flight- 
tested on the nose of a McDonnel -Douglas F-15 aircraft. A summary of the 
resulting noise and transition data has been- reported by Dougherty and Fisher 
(7). In this concluding report, Dougherty and Fisher found, for the range of 
Cp rms observed, that the data for transition Reynolds number, based on the 
product of local unit Reynolds number and distance from the nose to end -of 
transition , x^ 1- , appear to correlate with C p rms by the following equation. 

Re T ~< C P,™ S I” 25 <«> 

This relation, with the value of the proportionality constant suggested 
by Whitfield and Dougherty (24), is compared in Fig. 4 with some transition 
data obtained with the AEDC Cone in seven different tunnels (Dougherty and 
Steinle (1) and Mabey (25)) and a flight test at M - 0.80. 

oo 

The Dougherty-Fisher correlation indicates that the end-of-transition 
location, x T , is decreased by either increasing the tunnel noise or increasing 
the freestream unit Reynolds number. In order words, the effects of noise and 
Re m on x T are equivalent. However, their effects on measurements of C f or P p 
are not equivalent. Becker and Brown (26) have discussed the effects of 
turbulence on time-averaged pressures measured with Pitot probes. Since 




Re-r X 10“ 






turbulence causes fluctuations In the direction of the flow with respect to 
the probe's axis, the time-averaged pressure is reduced below the true total 
pressure. 1. Similarly, we have found that the tunnel noise, in the case of 
laminar boundary layers, also causes P p fluctuations and reduces P p 
measurements. This is equivalent to decreasing Re m . It- is important to 
distinguish between the effects of noise on_x T (which is the purpose of 
Dougherty, Steinle and Fisher's work) and noise effects on theoretical C f , or 
measured P p (which is the purpose- of this world. The two effects are actually 
opposite, Fig. 5. 


t 


This effect decreases as a wall is approached, since turbulence Is damped 
at an Impermeable wall. 




Figure 5. Effect of Tunnel Noise on Preston-Tube 
Measurements and Transition Onset 


CHAPTER III 
ANALYSIS PROCEDURE 
Experimental Data 

Although the AEDC Cone was tested In twenty- three different tunnels,, only 
the data, from the NASA Ames 11-Ft Transonic Wind Tunnel (27), 11-TWT for 
brevity, was considered for analysis In this study to develop subsonic wind- 
tunnel correlations. Table I lists nineteen subsonic wind-tunnel tests, and 
Table II lists nine subsonic flight tests which were selected for analysis in 
this study. The criteria for choosing a case for analysis are: 

1. The Preston-tube survey covers all three portions of the boundary 
layer. 

2. The flow angles a, e are very small. t. 

3. The freestream Mach number is less than unity. 

The distribution of static pressure along the surface of the sharp cone 
was measured only in flight. For wind-tunnel analysis, this pressure 
distribution is assumed to be defined by. the inviiscld-flow theory of Wu and 
Lock (20). Wu and Lock's predictions for the pressure coefficient along the 
surface of a 10-degree cone are shown in Fig. 6 as_a function of freestream 
Mach number. Measurements. of pressure coefficients together with linear curve 
fits from two typical flights are shown In Figs. 7 and 8. With this 
information and the known freestream conditions, the flow conditions at the 
outer edge of the boundary layer can be calculated. (For details see Ref. 8). 


,V eCeSSa 7 Since the bounda ry-1*yer code used in this 
u f? AN " 5 i 28 ’ wa „ s £ ound t0 be insensitive to changes in a, e. 
*h S °K V ? u ? s a > 0*5 and/or & > Q.25 8 have been shown to affect 
the beginning of transition, x B . Notice that the values tabulated in 
rabie II for a and_$ have an experimental uncertainty of ~ + 0.25°. 






TABLE II 


FLIGHT TEST CASES 


FLIGHT NUMBER M 


Re m x 10" 6 q„ (kPa) 



Note: a and & are time-averaged during a traverse 
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Figure 6. Inviscid Pressure Distribution on a 10° Cone (Wu & 










Correlation of the Data 

The fully-1 anil nar and fully-turbulent boundary layer computations are 
done using a computer program developed by Crawford and Kays (28) which they 
labeled STAN-5. The resulting distribution of skin friction and boundary 
layer properties are then correlated with the Preston-tube pressures. 

The form of the correlation equation Is derived from Eqn. v (4) using the 
parameters of Patel (11) and Quarmby and Das (12) but allowing the effective 
center of the probe to vary, l.e., 

y* * A(x*) 2 + Bx* + C (23) 


where 


X* ■ lO9 10 (Up K eff /V w ) 2 


(24a) 


and 


y* - log 1Q (U t T eff /v w ) 2 . (24b) 

Up and Y eff are defined as the longitudinal velocity and the height at which 
the theoretical total pressure (calculated by STAN-5) Is equal to the measured 
Preston-tube pressure at a given location on the cone surface. The 
coefficients A, B, and C are determined by a least-squares curve fit of the 
data. The results are presented and discussed in the next chapter. Fig. 9 
outlines the steps followed In the data analysis to obtain Preston-tube 
correlations. 
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Effective Reynolds Number Derivation 


Given the flight correlation in the form 

y* - Aj(x*) 2 + B 1 x* + C l (25a) 

and the wind-tunnel correlation in the form 

y* = A 2 (x*) 2 + B 2 x* + C 2 (25b) 

it is desired to derive an expression for the freestream unit Reynolds number 
in the wind tunnel when all other properties and parameters are equated 

between the two Eqns. (25a) and (25b) and the skin-friction coefficient 
predicted by the f 1 i gh t correlation is used. In other words, substitute the 

wind tunnel data into the flight correlation, solve for in flight, then use 
this value of together with the same wind tunnel data, except Re m> to solve 
for Re m , which is therefore the effective wind-tunnel unit Reynolds number, 

Re m,eff* 

The following identity relates the freestream conditions and can be 

derived using simple algebra, Abu-Mostafa (29). 

1.5881 x 10" 3 M Re 

— — T t - T « - no-3 = 0 (26) 

Thus, if only and q w are to be equated between wind-tunnel and flight 
correlations, then T^ must be allowed to change. This means that T Q also will 
change. Since it is desired to equate the values of the local Mach number, 
M g , between the two correlations so that the static pressure may also be 

equated, T e must therefore be allowed to change, and hence IL, p , T' and 



U . It can be shown that all variables other than T , T , T , p , T', U and 
P 08 o e e p 

U e can be kept unchanged without fixing Re m . Notice that T. Is assumed 

constant, T q g, along the cone for a given wind-tunnel case, but equals a 

different constant, T . , for the flight case. 

0,1 

Now, by substituting the definitions of x* and y* to both equations 
(25a, b) and subtracting one from the other to eliminate. Cf, the following 
equation Is obtained^- 

A 2 1o 9 1Q T 2 + A 2 + B 2 ) 1o9 10 T 2 " t4 F * (A 1 " A 2 } . 

+ 2F (B 1 - B 2 ) + (D x - D 2 + A l log^ 0 TJ + (4 + B l )log 1Q Tj ] = 0, 

where. . F - 1 og 10 (M p ( T R) * 5 Y eff /v w ). (27) 


This is a quadratic equation that can be solved for log^d^), hence T^, the 
effective local reference temperature in the wind-tunnel. T g g follows from 
the definition of reference temperature by Sommer and Short (15): 


£ — ■ 0.55 + 0.035M* ♦ 0.45 J*~ . 
e,2 e T e,2 


Then T^ ^ can be evaluated using the isentroplc relation 

,2 


->,2 


I + 0.2 M 


'e,2 


1 + 0.2 M 1 


(28) 


And finally Re „ (* Re cc ) can be calculated using Eqn. (26). 

m,z m,eff 

This procedure Is graphically outlined in Fig. 10, and the results of its 


application are shown In the next chapter. 
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CHAPTER IV 


RESULTS AND DISCUSSIONS 
Laminar Wind-Tunnel Correlation 

This and the following sections are discussed in detail in an earlier 

report by the authors (8). We will summarize the main points here. 

1. Laminar K eff values in the tunnel, plotted vs. R , indicated that the Re m 
= 9.8 x 10 measurements of Pp were in error as they were also in cases 
70.726 and 72.748. Correction was made according to NASA's directions 
(30) using case number 21.318 as a reference. The correction method is 
explained in point 5 below. 

2. The laminar wind-tunnel correlation, based on the shifted data, was found 
to be 

y* = -0.0103 (x*) 2 + 0.6653 x* - 0.5946, 

5.7 < x* < 6.3, 

0.30 <, £ 0.95, and (29) 

9.8 x 10 6 < Re m < 16.4 x 10 6 . 

3. The r.m.s. scatter of T f about Eqn. (29) is only 0.98%. This very low 
scatter demonstrates the importance of including a variable K e f f in the 
correlation. Without it, TT- „ me was found to be 4.93%. 



4. Eqn. (29) Is not universal. Its empirical coefficients contain 
information peculiar to the 11 TWT. In fact, no Preston-tube correlation 
is universal unless it properly models the wind-tunnel envl ronmental 


effects. 

5- Procedure for correcting laminar wind-tunnel data: 

The first objective is to align case 58.633 (M^ = 0.7 , Re m = 9.8 x 10 ) 
with case 21.318 (M w = 0.7, Re m * 13.1 x 10®) which is considered the 
reference. Then shift all the cases whose Re m = 9.8 x 10® accordingly. 

The second objective is to align case 70.726 (M^ = 0.7, Re m = 13.1 x 10®), 
with case 21.318, then shift case 72.748 (M = 0.7, Re m = 13.1 x 10®) 
accordingly. 

a. Evaluate of case 58.633 as the average of all R^ values in this 
case. Denote it by R KQ . 

T,DO 

b. Extrapolate the data in case 21.318 up to R^ g g. Use a French 

curve or do a least-squares curve fit of the data for case 21.318. 

c. Evaluate K ef f at R^ as given by the extrapolated curve; denote it 
by K ef f 21 * Also read K e ff at R x 53 as - 9 iven by case 58.633 (the 
original value). Call this value 

d. Compute AK eff & g = K ef<r 2 i - K ef f gg. This is the incremental 

adjustment of K eff for the Re m * 9.8 x 10® cases. 

e. Find aP q 5g = corresponding total pressure adjustment (from 

theoretical STAN-5 profiles). Add this increment, algebraically, to 
all P p measured values in case 58.633. 

f. Find AP^s for other cases in the Re m * 9.8 x 10® group which 

correspond to the same value of AK eff; above and shift these cases 


9 - 


by the proper increment of total pressure. 

The procedure for shifting cases 70.726 and 72.748 is similar to 
steps a-f above. 


Laminar Flight Correlation 

1. Original values of were plotted vs. R t> see Ref. (8). The plot 

strongly suggests that the flight data have random errors. 


2. The authors tried several correction approaches. The approach we 
recommended In Ref. (8) is based- on the. assumption that. the data of Flight 
#349.1400 are correct. The correction is done with the aid of the 
following equation which is a curve fit of the laminar, shifted wind- 
tunnel data. 


v eff 


-Q.273R 9 n 

= 2.865e T + 0.655 (1 - tr) u,i/J , 


0.60 < M < 0.95. 

CO — 


(30) 


Eqn. (30) was used to develop incremental adjustments to the flight values 
of K ef |r. It is plotted in Fig. 11 and is superimRUsed on the laminar, 
shifted wind-tunnel data. The correction procedure is outlined in point 5 
below. 


It was noticed that the correction of flight data resulted in changes to 

individual K - R distributions, compare Figs, 12 and 13. Since 

3K eff 3f> n 3P n 3P n 

— is a function of ~ and it is expected to change with 

movement across a profile. Fig. 14 clearly illustrates this idea for an 

actual case. 


36 



ro.7 





SYMBOL 



















3. The laminar flight correlation based on the corrected data was found to be 

y* * 0.05981 (x*) 2 - 0.1777 x* + 1.928, 

(31) 

5.6 < x* < 6.7, 

0.66 1 1 0.94, and 

6.9 x 10 6 < Re m < 9.2 x 10 6 . 

— m — 

The r.m.s. error in is only 0.37%. 

4. Eqn. (31), to the best of our knowledge, is the first-correlation in the 
literature for Laminar Preston-tube measurements in flight. 

5. Procedure for correcting the flight data: 

a. Let K e ff pp (R t ,Mj) = the value of at for the flight case 
with = Mj. Similarly define K ef p W j (R^.M^). 

Let = set of all R^ values in the flight case with = M^. 

Let {M^, M^}pp = set of all R^ values common between the two flight cases 
whose MJs are Mj and M 2 , i .e. (Mj ,M g } po = {M^p A {M 2 ) fd * 

Let R^M^, ^ 2 )pp = the average of all R^ values in {M^, M 2 )pp« 

Let t pg(^ 1*^2^ = ^eff ,FD^t*^ 1^ ” ^eff,FD ^t’^ 2^* ^ milarly define 

AK eff W j (M 1# M 2 ) ^ Refer now to Fig. 13. 

b. The reference case for all flight cases is flight #349.1400, i.e. M 2 * 
0.75. To shift a flight case {M^pp. first determine , 0.75} F p. 

If (Mj, 0.75}pp » 4 >, i.e., no R t values are shared by the two cases then 
we have one of two situations. 

{0.75}pp > {Mjlpp, in which case set R t (M^, 0.75)pp to be equal to 
the largest R^in {M.Jpp. 

(M^lpp > {0.75}, in which case set R^M^, 0,75)pg to be equal to the 
smallest R T 1n {Mj}pp. An example of such a situation is {0.66}pp, 
see Fig. 12. 
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Then, go to step_d below. 

c. This Is the case where {M^ 0.75 }p D Is defined $) such as {0.74}. So, 
calculate R t (M^, 0.75)^. 

d. F^nd K et£ ^ l) (R T ,M 1 ) and K effjFD (R T ,0.75) hence AK eff>FQ (Mj ,0.75) . 

e. Find K^ FF WT^x’^1^ and ^eff WT^t* hence yj(Mj,0.75) from 

curve-fit equation of K fiff vs. R t In the wind tunnel, such as Eqn. (30). 
Notice that R T (M X , 0.75) WT - 0..75) FD Also, AK eff >W r (M l *°* 75) wi11 

be negative if > 0.75. 

f— Calculate AK eff .shlft^llFD = ' incremental adjustment of K eff values In the 
flight Mj = AK eff(FD {M lf 0.75) - **< ef f s WT (Mj * 0.75). 

g. From the theoretical P Q profiles for the flight case M^, obtain 

^o.shlft^^FD which corres P° nds t0 ^eff.shlft^lVo at the Nation 
where R t = R t - This is the Incremental pressure adjustment for flight case 

Mi. 

h. For all points in {M x } FD , obtain P Pi# h1ft tM l l FD = p p< M l)FD ’ 

^o.shift^lVD* p p^ M l^FD be ^ n 9 the original, measured value of Preston- 
tube pressure. 


Laminar Effective Reynolds Number 

Based on Eqns. (29) and (31), the effective freestream unit Reynolds 
number was computed and plotted versus M m . The plot, Fig. 15 resembles the 
curve for noise data on the AEDC cone in the 11 TWT*, Fig. 16, and has a peak 
at M w ■ 0.70 - as does the noise. Actually, AR gff s (R© m eff - Re |n )/Re m 
correlated with noise by the following equation. 


t These data Include installation effects in addition to wall -generated 
noise. 
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0.07 t 


( 32 ) 


Aft 


eff 


6.2b (C 


P,rms) 


This supports the thesds that environmental effects In a wind tunnel can 
be calibrated by a single number, i.e., So, in order to measure the 

same average, theoretical skin friction coefficient, or the same measurement 
of Pp, in flight as is measured in the tunnel, the flight value of freestream 
unit Reynolds number should be increased to Re m eff . This effective Reynolds 
number will not necessarily equate the measurable values of C f . Indeed, the 
effects of noise on directly measured skin friction, if any, are unknown. 

The Transition Region 

Recall Dhawan and Narasimha's (21) intermittency function for 
transitional flow: 

-A5 2 (x) 

y(x) * 1 - e . (19) 


In order to be able to use Eqn. (19), x needs to be known for each case. 
Since measurements of y(x) are not available for this study, Eqn. (21) cannot 
be used. Another method was developed to calculate x as will be shown now. 

Calculation of A 

This method makes use of the available Preston-tube data__ Since it is 
assumed that the distribution of follows Preston-tube measurements (see 
Eqn. 4), one can assume that the location x T where P p peaks is the same 
location where C f peaks. Within the transition zone, the C f distribution is 
calculated using the y-function in the following manner: 


t The accuracy of this correlation is not very good since it does not 
include other environmental effects such as freestream turbulence 
intensity. Re m ff calibrates all these effects and not only noise. It 
should be noted, however, that noise effects are dominant in the ll-Ft 
Transonic Wind Tunnel, (1). 
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c f = (1 - y) C f>A + Y C f T , 


(33) 


where C f ^(x) is the local laminar skin friction coefficient if It were to 
occur at the given location x, and C f T (x) Is the local turbulent skin 
friction coefficient If It were to occur at x. The origin of the turbulent 
boundary layer is determined from the fully-developed turbulent flow at or - 
downstream from x E , the end-of-transltion location, as will be explained 
later. The value of x E corresponds to K = 4.0 (or y = 0.9986) as recommended 
by Dhawan and Narasimha (21). 

Differentiating Eqn. (33) with respect to x and evaluating at x-p yields 
the. f ol 1 owl ng rel atl on : 


dC 

dx 


-I = o = [(c c ) dY + ( dCf * T dCf **h + dCf **i 

l@x T 0 [(C f,T C f , £ ^x + c(x T~ )y + 3S ] x T 


(34)' 


A following formula for calculating Cf j Is reported by White (31) to be 
reasonably accurate. 


'f,T 


0.455 

S 2 Ln 2 [ Re 


w 


( V s - )- 5 3 

w 


Eqn. (34) Is also valid at x t , location of minimum P n . Solving for x R 
which appears In the definition of y. It was found thar x H « x. . 

Therefore, the value of x t is used from here on to designate the 
transition onset location. 
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Using Sommer and Short's model for S, a compress ibi lity factor (see Ref. 31), 
to correct for variable properties and Tetervin's (32) correction for 

/p /p' 

axl symmetric flow and making the approximation that — ~ -V the 

following equation can be derived. 




0.455 

c 37.8 v‘ J 


(35) 



J 


Here x y * distance along cone surface measured from the virtual origin of the 
turbulent boundary layer. It can be written in the form 

i 

i 

! 

\ - x - ax, (36) 

* 

r 

where Ax is the location of the virtual origin (see Fig. 17). It is now clear j 

j 

that Eqn. (34) can be solved for X if ax is known. The following section 

explains how this is done. I 

Calculation of Ax 

Eqn. (35) can be rewritten in the form ’ 




x v = 


- « - ^- 7 --- 8 - v - exp [ 
e 


0^455 P ’ 
C f,T ] 


0.5 


(37) 


So, all that is needed to calculate Ax is a reference C f T in the fully- 

developed turbulent flow at a location x - > x r . 

ref — E 

Crawford and Kays (28), who developed the STAN-5 program, state that 

their program's calculation of turbulent agreed with extensive measurements 

done at Stanford University. They used the following equation to effect 
gradual transition. 
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( 38 ) 


. , + „ 8e (x) - Re.(x R ) 9 

A + ’ A + ( ♦ (300 - A;)[l - sin ( 


Here A (x) is an effective sublayer thickness used in the van Driest damping 
model 


D = 


1 -Y + /A + 


(39) 


Fig. 18 shows a pjot of Eqn. (38) for a typical wind tunnel case. The damping 
coefficient is used in the Prandtl mixing length model for turbulent boundary 
layer calculations near the wall as follows.* 


t ■ k Y D, < = 0.41 


(40) 


And 


p w 



% U T 


dx 


Now, in Eqn. (38), it is assumed that 


(41) 


Re 0 (x E ) = 2 Re Q (x B ) . 

This was not found to be true at values of x £ * x g + 4X (recommended by 
Dhawan and Narasimha (21)). In addition, this transition model does not 


* The mixing length model is also the one used in this study to calculate 
the fully-developed turbulent boundary layer. 
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T » 1 

A + t = 300 
X F = Xref 




*■' TO? ■ llv-TT 


produce a peak in Cp at Xy. Instead, we used the followiog slightly different 
equation: 


A + - A*, ♦ C500 - A + t ][l - Sln( f • ( 4Z ) 

where Rep is the local length Reynolds number at a location Xp which is 

changed, so that a_peak in Cp occurs at Xy. This trial-and-error procedure is 

illustrated in Fig. 19. It is important here to mention that Eqn. (42) is not 
used as a transition model- Its sole role is to effect gradual transition so 
that the turbulent flow downstream is accurately computed. Indeed, when 

either of Eqns. (38) or (42) was used to simulate transition, the computed 

skin friction was found to be greatly underestimated as compared to the 
Dhawan-Narasimha model . 

To sum up, Eqn. (42) is used to prepare to compute turbulent flow, and 
hence obtain a good estimate of a reference value for at Xp or 

downstream. The location x pef >_ x £ is estimated from the Preston-tube data 
traces as the location downstream from Xp where the P p measurements exhibit a 
slope characteristic of fully-developed turbulent flow (see Fig. 20). 
However, this estimate of x pe p need not be precise, as long as it is 
sufficiently downstream from Xp, 

Using C f -j^ ref at x rg p and substituting in Eqn. (37), Ax may be 
calculated. Hence, X can be calculated from Eqn. (34). Thus, the y-functlon 
is now fully defined, and the C f distribution can be computed using Eqn. (33). 

In the above argument, it is assumed that White's formula, Eqn. (35), 
accurately calculates Cp y and/or x y . The authors have found, by trial and 
error, that C f y at x pef computed by STAN-5 and using Ax obtained from Eqn, 
(35) was always different from Cpj^ re p. Best results were obtained when the 
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virtual origin coincided with the transition point, i.e., 
x v (0) * x t = Ax. This finding was also reported by Dhawan and Narasimha. 
(21) Based on this finding an improved procedure to calculate x and hence y 
Is described next. 

The following variation of Wbite's—equation Is used in place of Eqn. 

(35). 



where C is a constant that_has a different value for each case and can be 
directly evaluated from Eqn. (43) at x pe1 r. 

Eqn. (43), then, together with 1ts_uieri vati ve w.r.t. x, the laminar STAN- 
5 calculations of ^ and its derivative w.r.t. x are substituted in Eqn. 
(34) to solve for X and hence y. 

The Transition Correlations 

In order— -to completely define the correlation parameters x* and y*, 
theoretical velocity and total pressure profiles in transition need to be 
computed to obtain Up(x) and K g ^^(x). These profiles may be calculated using 
the y-function in a manner similar to skin friction, Eqn. (33). 

U(Y) = (1 - y) U a (Y) v Y U t (Y), (44) 

T(Y) - (1 - y) T a (Y) + yT t (Y). (45) 

From these two profiles, calculate P 0 (Y) as follows: 



M(Y) « U ( Y ) /( 555 /RyT) 


P q (Y) = P W C1 + 0.2 M 2 (Y)] 3 * 5 


( 46 ) 


Initial- profiles for turbulent flow computation can _be_ obtained by 
rescaling available fully-developed turbulent profiles (at x pef ) using edge 
velocity and boundary-1 ayer thickness at the initial location which can be 
estimated using Musker's equation, Musker (33), as follows: 


0.41 U 


fi = e 


3.0504 


'' U X 


at x 


initial * 


(47) 


where U e = ( ^ ) °* 5 at x in1t1a1 . C fT at x init1al can be calculated 
using Eqn. (35) with x v = Xinitial - ax and all properties evaluated at 
x initial downstream from x^.. 

Values of Up and can then be computed by interpolation of measured 
Preston-tube pressures in velocity and total pressure profiles given by Eqns. 
(44 and 46). 

Based on the above analysis the transition correlations for the original 
data are: 


Wind Tunnel : 

y* = 0.06935 (x*) 2 + 0.02795 x* + 0.9678, 

5.2 < x* < 6.3 , 9.8 x 10 6 £ Re^ £ 16.4 x IQ 6 , 0.30 £ £ 0.95, 

* = 2 - m » and (48) 


(49) 


y* « 0.02094 (x*r + 0.5988X* - 0.7112, 

5.5 < x* < 7.1 , 6.9 X 10 6 < 9.2 x IQ 6 , 0.66 < M < 0.94, 

C f,rms * 3 - 64 «* 

A plot of. Eqn. (48) with the superimposed wind-tunnel data appears in n q . 
21. Fig. 22 Is a plot of C f scatter about Eqn. (48). Figures 23 and 24 
illustrate the same for the flight data. 

Not all — the available data in transition are included in the above 
correlations; only the points at which x* and_y^ are proportional are included 
(These amount to slightly more than 60% of the total number of points in the 
transition region.). This requirement is suggested by. the basic Eqn. (4). 

Figures 25 and 26 are plots of -transitional values of K gff vs. in the 
wind tunnel and flight, respectively. Notice that the data, again, indicate 
large errors in the flight tests. Before discussing how these errors are 
corrected, we first present the results from the analysis of turbulent data. 

The effective Reynolds number distribution based on Eqns. (48) and (49) 
is shown in Fig. 27. It does not correlate with noise. This situation may 
change after correcting the experimental data. 

The Turbulent Region 

The procedure, which is described two sections earlier, for estimation of 
a reference C f>T provides an accurate and complete method for theoretical 
computations of C f , velocity and enthalpy profiles in the turbulent flow 
region. Therefore, all the information needed to define x* and y* for this 
region is available. 
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Figure 24. Scatter of Transiti 
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Figure 25. Distribution of Effective Probe Height as Determined from the 
Original Transitional Wind Tunnel Data. 





Figure 26- Distribution of Effective Probe Height as Determined from the Original 
Transitional Flight Data. 
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Figure 27. Distribution of Transitional Effective Reynolds Number Based on Original Data 







The wind tunnel data are corrected In a manner similar to the laminar 
data, viz., by referencing all cases to case 21.318 (M, * 0.7, Re,, - 13.1 x 

10 ’ q * " 26-3 Un,1ke the laminar data, the Re„ = 9.8 x lo 6 cases 

already form continuous curves of K eff versus R x , Fig. 28 . so, the only cases 
which are shifted are cases 70.726 (M„ - 0.7, Re„ . 13.1 x lo 6 , q . 25 8 
kPa) and 72.748 («. » 0.8, Re„ - 13.1 x 10 6 , - 29.0 kP„. 

Similarly, the flight data are corrected in the same manner as the 
laminar data, see outline at the beginning of this chaster. 

The turbulent correlations without corrections are found to be: 

Wind Tunnel : 

y* = 0.02337 (x*) 2 + 0.5Zl5x* - 0.6202, 

5.1 < x* < 6.9, 9.8 x 10 6 < Re m < 16.4 x IQ 6 , 0.30 < M w < 0.95, (50) 

C f,rms = l ' m » and 
Flight : 

y* = 0.007512(x*) 2 + 0.7749x* - 1.272, 

6.0 < x* < 7.7, 6.9 x 10 6 < Re„ < 9.2 x 10 6 , 0.66 < < 0.94, 

6 f,rms = 1*10S • 


Eqns. (50, 51) are plotted with the data In Figs. 29 and 30. The scatter of 
C f is shown in Figs. 31 and 32. Figures 28 and 33 show the distribution of 
K eff vs - r t - Notice that the relative positions of different flights in 
Tig. 33 Is the same as shown In Fig. 21 of Ref. (8) which Is reproduced In 

Tig. 12. This suggests that the same correction procedure can be successfully 
applied. It was Indeed as will be shown shortly. 

The effective Reynolds number distribution based on Eqns. (50) and (51) 
is Shown in Fig. 34. Again, it does not look like the noise curve. Fig. 15, 
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Figure 28. Distribution of Effective Probe Height as Determined from 
the Original Turbulent Wind Tunnel Data. 
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Figure 31. Scatter of Turbulent Skin Friction Coefficient About 
Correlation for Original Wind Tunnel Datk. 
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gure 32. Scatter of Turbulent Skin Frictipn Coeffi 
About Correlation for Original Flight 
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riginal Turbulent Flight Data 








which may be caused by the errors In the experimental data. The correction 

procedure used to correct the laminar data should result In a aH „ 

eff 

distribution which Is closer to the noise distribution, as can be observed In 
Fig. 36 of Ref. (8) reproduced In Fig. 15. 


Results After Data Corrections 

The Turbulent Region: 

The turbulent wind tunnel data after shifting a subset of It as explained 
before are shown in Fig. 35. The correlation is given by 


y* = Q.02282(x*r + G.5782X* - 0.6409, 

5.1 < x* < 6.9, 

9.8 x iO 6 < Re < 16.4 x 10 6 , 

— m — 

0.30 < M < 0.95 and 

*— CO — 

C f,r m s ■ 1 - 20 * • 


(52) 


Notice that there is no significant change to the correlation coefficients and 
accuracy since the shifting was minor. Eqn. (52) is shown in Fi_g. 36 with the 
data and the scatter of these data about Eqn. (52) is shown in Fig. 37. 

The corrected flight data appear in Fig. 38. Notice, again, that the 
distributions of K e f f versus R t for Individual cases has been altered by the 
correction procedure. The flight correlation is given by 

y* = O.GQ5586(x*) 2 + 0.7723 x* - 1.1867, 

5.45 < x* < 6.30, 

6.9 x 10 6 < Re m <. 9.2 x 10 6 (53) 

0.66 < M < 0.94 and 

— » — 

f , rms *■ 0.65% . 


72 


<i 



Distribution of Effective Probe Height as Determined 
1 from the Shifted Turbulent Wind Tunnel Data. 









Turbulent Correlation for Shifted Wind Tunnel 






Scatter of Turbulent Skin Friction About Correlation for 
Shifted Wind Tunnel Data. 






SYMBOL 









This equation and the corrected data are shown In Fig. 39 and the data scatter 
are shown In Fig. 40. 

Based on Eqns. _(52) and (53). the AR fiff distribution Is shown in Fig. 
41. The distribution does-not bear any resemblance to noise characteristics , 
Fig. 16. This -means that despite the data correction, the information 
contained in them and their correlations are not sufficient to extract the 
expected Re m>eff information. The reason for this is the added complexity 
that is not present in the laminar analysis, namely the vortlcity fluctuations 
m the boundary layer. As other investigators have found, these fluctuations 
are large enough to dominate pressure fluctuations caused by background noise 
and thus eliminate their effect on Preston-tube measurements. For example, 
Whitfield and Dougherty (24) reported the effect; of background noise on 
transitional and turbulent boundary layers on the AEDC cone in four transonic 
wind tunnels. They noted that each of these tunnels had an acoustic resonance 
near M w = 0.8, but that "the frequency components coming into resonance in 
these slotted -wal 1 tunnels were so low (< approximately _2Q0 Hz) that the cone 
boundary layer was Jnsensjtlye to them and their influence on transition was 
nil." Weeks and Hodges (34) also concluded that even at noise levels up to 
C p,rms = 8% U1t was not possible to identify any effect of the noise itself on 
the boundary layer, and it is concluded that the acoustic disturbances 
generally found in the working sections of transonic wind tunnels are unlikely 
to exert measurable Influence on the development of turbulent boundary layers 
on wind-tunnel models - at least for mild pressure gradients." Furthermore, 
Raghunathan, et al. (35) showed that turbulent skin friction coefficient was 
haraly affected by noise levels up to C p>rms = 2%. Based on these findings, 
the value of AR eff for turbulent data is expected to be zero for flight and 
wind tunnel cases with identical freestream flow conditions. Wind tunnel case 
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Figure 39. Turbulent Correlation for Corrected Flight Data 











Figure 41. Distribution of Turbulent Effective Reynolds Number Based on Corrected Data 







for 


#56.631 and flight case #333.1354 have similar flow conditions, and aR 

ef f 

these conditions Is indeed near zero, see Fig. 41 at M = 0.90. 

00 

As noted before, Becker and Brown (26) showed that pressure fluctuations 
decrease the measured Preston-tube pressure. Pressure fluctuations may be 
caused by background noise and/or by Internal boundary layer turbulence. 
Since voctlcity fluctuations in a laminar boundary layer are negligible, 
background noise and freestream turbulence dominate in this region, and the 
data analysis described in this report permits the calibration of these 
environmental effects. In transitional and turbulent boundary layers, on the 
other hand, internal fluctuations are dominant and background noise has no. 
effect on the measurement of P p and, therefore, cannot be calibrated. The Re 

g 

s 9.8 x— 10 data show the greatest deviation of aR eff from zero. Fig. 41. The 
reason is that this group of data has the greatest experimental uncertainty in 
Pp measurements, and it is the reference for correcting the flight data 
(i .e. , values at which- correction is made correspond to a wind tunnel 
Reynolds number of 9.8 x 10®. Thus, the corrections for the low Reynolds 
number, turbulent flight data appear to have been inadequate. 

The Transitional Region; 

In order to Insure the continuity of the K fiff distribution during transition, 
the ^q, shift 1ncrements us ed in the correction of flight data must vary 
gradually from the AP values used in the laminar correction and those 

used in the turbulent correction. We used a linear variation in the following 
form: 


^o.shift 


c AD 

o, shift,* 


x - x t 

xfr^ (Ap o, shift, T 


AP o, shift,*)* 
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F19. « shows the continuous K #ff - R t distribution for esse 19.289 In the 
three regions of the boundary layer. The results after shifting the wind 
tunnel datajnd correcting the flight data are as follows. 

Wind Tunnel: 

y* = 0.7814(x*) 2 - Q.07967x*_±_1.2936, 

5.25 < x* < 6.30 , 

9.8 x 10 6 < Re m < 16.4 x 10 6 , 

0.30 < M < 0.95 (MJ ... 

— » <30 — * 

C f,rms = and 


Flight: 


y* = 0.09131 (x*) 2 - 0.2596x* + 1.9066, 
< x* < 

6.9 x 10 6 _< Re m _< 9.2 x 10 6 , 

0.66 < _< 0.94 and 


'f.rms 


2.31S 


(55) 


Eqn. (54) with the wind tunnel data and. their scatter are shown In Figs. 43, 

44. Fig. 45 shows the K fiff - R t distribution. Figs. 46 through 48 show the 

same for the flight data. 

Based on Eqn. (54) and (55), the AR eff distribution Is shown In Fig. 

49. As expected the distribution cannot be correlated with noise effects for 
the sane reason discussed In the turbulent analysis last section. 

Furthermore, the authors (8) have shown that the extent of transition, x - 
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Y* = 0.7814(X*r -0.07967X^ + 1.29 
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Figure 44 Scatter of Transitional Skin Friction Coefficient About 
Correlation for Shifted Wind Tunnel Data. 




stribution of Effective Probe Height as Determined from the Shifted 
Transitional Wind Tunnel Data. 
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Figure 47. Scatter of Transitional Skin friction Coefficient About Correlation 
for Corrected Flight Data. i 




Figure 48. Distribution of Effective Probe Height as Determined from the Corrected 
Transitional Flight Data ! 
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x t , is larger in flight than in wind tunnel tests with the same flow 
conditions. This means that the transition process requires a larger distance 
in flight than in a wind tunnel, and the rms values of P p indicate the laminar 
break-down in flight is more violent and, hence, creates larger vorticity. 
Indeed the rms fluctuating Preston-tube pressure coefficient in flight 
#333.1354 is nearly twice that In wind tunnel case #56.631. (These are the 
two cases with similar freestream conditions). 

We conclude, therefore, that the calibration of wind tunnel envi ronmentaJ 
effects on Preston -tube_ measurements or theoretical skin friction by an 
effective freestream unit Reynolds number can only be achieved by analyzing 
the laminar data as described in this report. 
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CHAPTER V 

SUMMARY AND CONCLUSIONS 


A new procedure has been developed which uses Preston-tube data from wind 
tunnel and flight tests_a f __the AEDC Transition Cone to compute an effective 
unit— Reynolds number for transonic wind tunnels. The resulting effective 
Reynolds numbers are based on the requirement that the average Preston-tube 
pressure for a given type of boundary layer be equal in the wind tunnel and 
flight for a given and q w but differing Re m . The procedure has been 
applied to laminar, transitional, and fully-developed turbulent boundary 
layers by using wind tunnel data obtained in the 11-TWT. The results for 

laminar boundary layers indicate that noise in the 11-TWT causes Preston-tube 
pressures to be low compared to the values that exist in flight for the same 

Re m> and 9«* Thls results in the effective unit Reynolds number being 
— 1 - 9 h ? - r . than the reference or operating value by approximately 6.5%. Thus, in 
order to increase the laminar Preston-tube pressures, obtained in the 11-TWT, 
to match the corresponding flight data, it is necessary to increase the tunnel 
unit Reynolds number by 6.5%. 

This unit Reynolds number trend is opposite to what is found in the 
technical literature on the effects of noise on boundary layer transition. In 
that context, transonic wind tunnel noise is known to promote early transition 
and is frequently viewed as being analogous to an increase in unit Reynolds 
number. With this perspective of matching the location of transition, 
transonic tunnels are though to have "effective" Reynolds numbers that are 
somewhat higher than the operating value selected by the tunnel operators. 
However, if for example a transonic tunnel is operated at a lower unit 
Reynolds to achieve matching of flight values of transition location on the 
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AEDC cone, one would not expect a match in drag values. In fact, the lower 
tunnel Reynolds number would result in lower skin friction within both the 
laminar* and turbulent portions of the boundary layer. 

Unfortunately, actual measurements of skin friction were not performed in 
either. -the wind tunnel or flight tests. Thus, the authors are unable to reach 
any definitive conclusions as to the effects of tunnel noise on the actual 
laminar skirufriction per se. 

Ihe-basic achievements of this study are summarized below. 

1. The law-of-the-wal l_is a valid w av to correlate Preston-tube data in the 
form of Eqn. (4) or Eqn. (23). 

2. The effective height of a Preston tube is not fixed.. It varies with 
u T h/v w , aspect ratio and the position of the probe with respect to 
the wall. Chapter II, p. 13. 

3. Including a variable K fiff in the correlation substantially improves its 

accuracy. Chapter IV, p. 34. 

4. — Blotting vs. R^ permits the detection of errors in the experimental 

data, .Chapter IV, p. 34. 

5. The wind-tunnel data can be used to correct errors in. P p measurements in 
the flight tests. A systematic correction procedure was developed and 
successfully applied to the flight data, Chapter IV, p. 41. 

6. Preston-tube correlations for laminar, transitional and turbulent data 
were obtained both for the wind-tunnel and free-flight tests. The flight 
correlations, Eqns. (31), (53) and (55), are the first of their kind in 


t 


This assumes that noise does not change the steady-state laminar skin 
friction in any significant amount. 


the literature. 

7. A semi -empirical method has been developed to define and calculate an 

effective Reynolds number which calibrates environmental effects on 
Preston-tube measurements in wind tunnels, Chapter III, p. 31. 

8. A computational model for the transition zone can be devised using 

laminar and ful ly-devdloped turbulent calculations (or measurements) of 

and transitional Preston-tube pressure measurements without the need 
for hot-wire measurements of the intermittency factor, y, Chapter IV, p. 
44. 

9. The. virtual origin of the turbulent boundary layer on the AEDC cone 

coincides with the onset of transition which is found to occur at the 
location of minimum P p , viz., x t . Chapter IV, p. 46. 

10. Experimental Preston-tube pressure measurements appear to have smaller 

errors in the turbulent portion of the boundary layer than in the other 

two portions, compare Figs. 25 and 28. 

11. The_ effective freestream unit Reynolds number distribution obtained from 

the analysis of laminar data is apparently correlated with noise data on 

the AEDC cone, Eqn. (32). Therefore, calibration of environmental 

effects in a wind tunnel can be done by calculating using laminar 

measurements of Preston-tube pressure. Best results are obtained when 

the freestream flow conditions, M , Re , and q , are the same in the 

eo m ® 

tunnel and in flight. 

12. The analysis of transitional and turbulent Preston-tube data may not be 
used to calculate AR gfrf since vorticity fluctuations in the boundary 
layer make it insensitive to background noise. The derived AR^'s from 
these data do not calibrate the tunnel's flow quality, but rather reflect 
the effect of internal vorticity fluctuations on P p measurements, Chapter 
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Vortlclty and pressure fluctuations In transitional boundary layer flow 

ace -l arger in flight than in the 11 -TWT for similar freestream 
conditions. Chapter IV, p. 91 . 

A traversing Preston tube is insufficient, by itself, to calibrate the 
effects of transonic wind-tunnel noise on skin-friction. The Preston- 
tube data must be supplemented with direct measurements of skin friction 
if this objective is to be achieved. 



CHAPTER VI 

RECOMMENDATIONS - 

The calibrating procedure described in this report may be used to 
calibrate environments in other transonic wind tunnels, especially those 
tunnels, where the AEDC cone was tested. 

The author recommends that skin friction be measured directly and used in 
conjunction with Preston-tubes In future wind tunnel_and flight tests. This 
will permit the described calibration procedure to reveal the effects ^of. 
noise. If any, on skin friction drag. 

Care should be taken in measuring Preston-tube- pressure in future 
experiments. Every effort to prevent probe twisting and lifting will reduce 
experimental errors especially in the flight tests. The gain factor and the 
reference pressure for the transducer should be accurately recorded. 
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